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We describe an opto-mechanical system in which the coupling between optical and mechanical de- 
grees of freedom takes the form of a fully quantised third-order parametric interaction. Two physical 
realisations are proposed: a harmonically trapped atom in a standing wave and the 'membrane in 
the middle' model. The dominant resonant interaction corresponds to a stimulated Raman process 
in which two phonons are converted into a single cavity photon. We show that this system can ex- 
hibit a stable limit cycle in which energy is periodically exchanged between optical and mechanical 
degrees of freedom. This is equivalently described as a parametric self-pulsing. 

O 

^ . INTRODUCTION 

^ ' The ability to couple optical and mechanical degrees of freedom has driven innovation in many areas of modern 
physics, from the earliest demonstrations of laser cooling of atomic motion[l] to the pioneering theory of Tombesi 
and CO workers[^ to recent experiments using macroscopic mirrors coupled cavity fields[3]. In the latter example, 
^-H , most investigations are based on the radiation pressure coupling between the photon number in a cavity and the 
' displacement of a mechanical resonator, typically a mirror, which is linear in that displacement Q ■ 
, In this paper, by contrast, we describe a more complex dynamical situation in which the photon number in the 
O i' cavity couples to the squared displacement of a mechanical resonator. We will consider two physical implementations, 
^Jj , one based on harmonically trapped atoms interacting with a standing wave in a cavity via the off-resonant dipole 
^ interactionfH', ^ , and the 'membrane in the middle' approachj?] . This means that the cavity field experiences a 
^ phase modulation at twice the mechanical vibrational frequency. If this modulation frequency is larger than the 
^H' cavity line-width, we can resolve these sidebands spectroscopically and drive the cavity on the first red sideband. 
In that case we can linearise the interaction with respect to the cavity field (not the mechanical system) and the 
dominant resonant terms then take the form of a fully quantised third-order parametric interaction. We show that 
this system can exhibit a stable limit cycle in which energy is periodically exchanged between optical and mechanical 
' degrees of freedom. A semiclassical analysis, based on the centre manifold, enables us to calculate the amplitude and 
period of these oscillations. A first approach to a full quantum treatment enables us to calculate the phase diffusion 
' constant on the limit cycle. 

, Domokos and Ritschf&l describe a system in which one or more atoms are trapped by the optical dipole force of a 
• standing wave in an optical cavity and, as they move, modulate the standing wave amplitude. Zippilli and MorigiQ 
, consider a variation on this scheme in which the atoms are independently trapped by a harmonic potential provided 
by another focussed optical dipole trap or an ion trap. We consider a similar model but, like Domokos and Ritsch, 
include the back-action of atoms on the standing wave of the optical cavity. If the atoms are trapped at a node of the 
cavity field, the Hamiltonian for the system is 



> 



0^ 
O 



X 



H = hujca^a + hvb^b + ;i(e*ae*"* -|- £^0^6"*"*) -f hGa'^ a{h -f b^f (1) 

where lOc is the cavity optical frequency, v is the mechanical trapping frequency, tc represents a coherent cavity driving 
field at frequency u, and the coupling constant is given by 

G = ^ (2) 

where g is the single photon Rabi frequency for the atomic dipole transition, A is the detuning between the atomic 
dipole resonance and the cavity field. The Lamb-Dicke parameter is 
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and k is the wave number for the cavity field. Typically the atom is harmonically bound with a typical trapping 
frequency in the range v ^ 10* — 10^ Hz. This is much smaller than the optical frequency, but can be of the order 
of the detuning between the cavity field and the cavity driving field. In an interaction picture at frequency co the 
Hamiltonian may be written as 

H = hSa^a + hvb'^b + n(e*a + e^a)) + hGa)a{b + h'^f (4) 

where 5 = ujc — uj. 

In the absence of the opto- mechanical coupling (G = 0), the cavity field will reach a steady state which is a coherent 
state with amplitude {a)ss — ol with 

(5) 



K/2-i5 

where k is the photon decay rate for the cavity field. In the presence of weak opto-mechanical coupling we proceed 
by transforming to the displacement picture to remove the coherent state component of the optical mode, 

H ^ D\a)HD{a) (6) 

Then 

H = h5a)a + hvb'^b + hG\a\'^{b + b'^f + hG{aa) + a*a){b + b'^f (7) 

We will further assume that the mechanical frequency is larger than the cavity line-width, > k, so that the optical 
cavity can be driven on resolved sidebands. We then transform to an interaction picture via the unitary transformation 
U = exp[— z(5ta^a — ivtb^b], and assume the parametric resonance condition 

5 = ±2i/ (8) 

Keeping only resonant terms we have two different kinds of resonant interactions 

Hr ^ hG{aah^ + a*ab^^) + 2hG\afb^b (9) 
when the laser driving the cavity is tuned to the red (lower frequency) of the cavity frequency uj = tOc — ^.v, and 

Hb = hG{aa^b^^ + Q*ab^) + 2hG\a\'^b'^b (10) 

when the laser driving the cavity is tuned to the blue (higher frequency) of the cavity frequency lo — uj(. + 2v. These 
Hamiltonians describe two-phonon Raman processes. In what follows we will absorb the Stark shift of the mechanical 
frequency (proportional to |ap) into the definition of the mechanical frequency. 

Thompson et al.[7] describe an opto-mechanical system comprising an optical cavity containing a thin dielectric 
membrane of SiN placed between rigid high-finesse mirrors. The Hamiltonian describing the interaction between the 
optical mode and the mechanical displacement of the membrane is given by [7| 

H = hujcix)a^a + hub^b + ft(e*ae"^* + ecfl^e"''^*) (11) 

where a, refer to the optical mode and b, b^ refer to the mechanical resonator and we again include a coherent 
driving field on the cavity and 



-ib + b^) (12) 
Zvm 

If the membrane is placed at a; = 0, an extremum of lUc, we approximate 

H = hwco^a + hvb^b + h{elae"^^ + £^0^6"'"*) -f hGa) a{b -f b^f (13) 

where uJc = ijJc(0) and 



G = 



n & 



|2, 



(14) 



x=0 

The Hamiltonian in Eq. (jl3p is the same form as the previous model in Eq. ([1]) . In a similar way we can now consider 
two resonance conditions depending on whether the cavity field is driven on the red side band or the blue side band. 
In this paper we will consider only the red sideband case given by Eq. ([9]). 

The possibility of realising two phonon Raman processes in this system has been noted by Thompson et alQ- 
For the parameters given in Thompson et al.Q, k ^ 10"* — 10^ s^^, 7 ~ lO^s^^ while the mechanical frequency is 
J/ ~ 8 X 10^s~^. Thus it should be possible to achieve the resolved side band regime for an optical cavity of sufficiently 
high finesse. 
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RESPONSE TO MECHANICAL DRIVING 



We consider a cavity mode coupled to a mechanical resonator via the interaction Hamiltonian in Eq. ([9]). This 
describes a situation in which two times the mechanical frequency, 2z/, is larger than the cavity decay rate, k, so that 
we can drive the cavity on the first red sideband. We now consider the response of the cavity field to a coherent driving 
force on the mechanical resonator. This is the analogue of the nonlinear optical process of second harmonic generation 
in which excitations in the driven mode at frequency v are converted to a single excitation in the un driven mode at 
frequency 2z^. In the case considered here however the transitions are Raman-like processes in which two phonons are 
up-converted into one cavity photon by absorbing cavity pump photons at the sideband frequency lOc ~ ^v. 

Moving to an interaction picture and neglecting counter rotating terms in the interaction between cavity field- 
mechanical resonator, we obtain the interaction picture Hamiltonian, 

Hi{t) = fix {a^b^ + a{h^ f) - + (15) 

where, for simplicity, we assume that we can choose the phase of the cavity sideband driving, ec, to make a real and 
we define 

X = Ga (16) 

Note that while G may be small Ga can be large. This is very similar to the strong coupling limit recently achieved 
for a linear opto-mechanical systemQ 

In addition to the explicit optical and mechanical degrees of freedom included in the Hamiltonian, we need to model 
dissipation and associated noise sources. The resulting master equation ( in the interaction picture) for the total state 



of optical plus mechanical oscillators is then[ll| 



^ = ~^[Hiit),p]+KV[a]p + jin+l)V[b]p + jnV[b^] (17) 

where k, 7 are the amplitude damping rates of the optical and mechanical oscillators respectively, n is the mean 
bosonic occupation numbers for the mechanical bath. The super operator V is defined by 

V[A]p = ApA^ ~ i {A^Ap + pA^A) (18) 

The master equation, Eq. (fT7)l , has been previously used to describe sub-second harmonic generation for two quan- 
tised fields interacting in a medium with a significant second order optical nonlinearitvflU . The first step to under- 
standing the behavior is to consider the semiclassical equations of motion. These exhibit a rich behavior including 
fixed points, and limit cycles. We will follow the treatment of Drummond et al. [l3|. We proceed by first analyzing 
the classical equations of motion for the system. 

Semiclassical dynamics 

The classical non-linear equations of motion for the red sideband interaction are, 

together with their complex conjugate counterparts. Here a and (3 are the complex amplitudes of the cavity field and 
the mechanical oscillator respectively. Setting x = 1, as it simply scales time, and rewriting in real form gives 

Pr = 2{(3^ar- f3rai)-^(3r, (20) 

A = 2(/3,a, -f Aa,) - |a - e (21) 

dr = ~2(3r^3^- '^ar, (22) 

= - - (23) 
This system has one critical point existing for all e. 
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The real parts of the variables at the critical point are zero; PrO = 0, a^o = 0. The imaginary part of /3 is then 
given by the solutions to the following cubic 







(24) 



Note that if e is assumed to be positive then /3jo is negative. (The other case (e < 0, (3io > 0) can be obtained from 

symmetry.) The imaginary part of a is given in terms of the solution on the cubic; ajo 
The linearized matrix about this critical point is 



/ -2aio - ^ 



V 





-2Ao 




2Ao 

2 



\ 



2a,o - ^ 2/3,0 





-2Ao 







■f / 



Now at a possible Hopf bifurcation two of the eigenvalues of this linearized matrix are pure imaginary. This can occur 
here when a,o ~ — '^'^'^ 



^ '-, which gives 

2 _ k(k + 7) 



€ = eh = - 



Ao(k + 27) ^ Vk(k + 7)(k + 27) 
2 aV2 



(25) 



Note that for k » 7 then eh ~ • 

We have assumed a situation where the Hopf bifurcation is analyzed as e is varied away from and the other 
parameters are fixed. So the calculation of the center manifold at the bifurcation is performed for e = Ch- The 
linearized matrix at the Hopf bifurcation is 



f 2Ao \ 

-7-1 Ao 

-2Ao -f 

-2Ao -f / 



where = Pmh = ~'\/^^^^- This has eigenvalues 



Aid 



±i 



. ^k{k + 27) 



A2± 



K + 7 _|_ ^ a/2k:(k + 7) -7^ 



The equations of motion after perturbing off the critical point at 
(0, Pioh, 0, a2io = 2^), are 



/ f3r\ 


/ 


dr 








V / 


V 



2 






2/3io/» 

K 

2 








-7-f 2Ao/. 



















+ 


/ 




V 



(26) 



The center eigenspace is the (/3r, Q^r) space and the center manifold is tangent to the {Pr, ctr) space. Say it is given, 
at least locally by /3j = /ii(/3r, cir) and a, = h2{Pr, cir), then to lowest order in Pr and the center manifold is given 
by 



Pi = hi{Pr, ar) ~ AiP'^ + BiPrar + Cia^ 

O'l = h2{Pr, Cir) ~ A2PI + B2Prar + ^20:^ 



(27) 
(28) 



Substituting back into the equations of motion gives values for Ai, Bi, Ci. 

2 - 2V2 v^<kT7)«(27k^ + 927k2 + 96^72 + 167^) 



D 



B, = 



4«;2(11k2 + 34K7 + 327^) (27 + 3/t) 
D 



5 



-(4(27 + Sk))V2^/k{k + j)k{k'^ + 2k-/ + 47^) 
= D 



2k{5k^ + 247«;2 + 32K-f^ + 167^) (27 + 3/t) 

I2 ~ 



Bo 



D 




-(8(27 + 3k))V2v'k(k-F 


-7)k2(2k + 57) 


D 




(8(k + 27))(5«; + 27)k(k 


+ 7) (27 + 3k) 


D 



(^2 

where D = (k^(47 + 3k)(327^ + 96k;7^ + 72«;^7 + 17k^)). Thus the equations of motion on the center manifold, for 

e = Ch, are 

2{-A2(3^ + (Ai - B2)a./32 + (Bi - C2)/3,a2 + C^a^) 



(29) 



where /3io/t ^ - y ^^^^^tli. 

Transforming to normal form 



2 2 



which has the inverse 



u = — , = H , = , (31) 

Vk(k + 27) 2VK(K + 27)/3^o;i 



The equations of motion in the normal form variables are 



where Ai, Bi, Ci, Di are functions of k and 7 and can be calculated from the equations above. 

These equations are valid for e = Ch, however they can be extended to include the case where e = Ch + Ae, by 
considering the variation in the trace of the linearized matrix with e. Recall that the linearized matrix for (/J^, ctr) at 
the critical point is 

-2aio - 2 2Ao 
~2/3jo — I 

2/3? 

where Um = —. The trace of this matrix is a function of e; 
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since must also be treated as a function of e via the cubic 



To a first approximation in Ae 



On calculation 



dtvCLCG 

trace{e) = trace{eh) ^ ^^|(£=eh)Ae 



dtrace y/8K{K + 7) 

de " 3k + 47 

which means that the linear system for e ~ eh in polar coordinates is 



r = dAer where d = ^ , ^ ~^ 

k(3k + 47) 

Finally the effect of the nonlinear terms is to include an term, whose coefficient is a function of the nonlinear 
coefficients given above(for details sec [9]) 

k2 (k + 7) {mn^ + 4907^3 + 808^272 + 512^73 + 1287-*) 
" ^ 4(128^274 + 480^373 + 51k6 + 284^57 + 576K472) 

Note that for k > 7 then a « 

Since d > and a < there is a supercritical Hopf bifurcation creating a stable limit cycle, which exists for A small 
and positive. In polar coordinates this gives the approximate system for the radial equation as 

A 3 • - 2 dAe 

f — dAer + ar which has non zero solution r — . 

a 



This gives the amplitude of the limit cycle: A = \ . The frequency of the resulting periodic orbit is 

a 



For K ^ 7 the approximate amplitude and frequency of the resulting oscillation are 



^ 1 /l36%/2Ae n , ^, 

A.-^^^ .h^-. (34) 

However caution should be used when considering this equation as it is the amplitude in the transformed variables. 
As we will see later on the amplitude in the original variables is multiplied by k. The solution on the limit cycle in 
normal form variables is 

u = Asin{u;iit + (35) 
V — Acos{LJht + (j)). (36) 



In the original variables this becomes 



and 



Pr = 2p,ohAcos{uht + (l3) (37) 
ar — u;/iv4.sin(cj/ii + 0) — cos{uJhi (1^) (38) 



Pi = Pio + hi{2PMhAcos{ujht + 4>), LOhAsiii^LUht + 4>) — cos{u!ht + (f)) 
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oil = Oio + h2{'2(3iQhAcos{ujht + (/)), ti;,iv4sin(cj/ii + 0) 2" cos(wftt + (/))) 

However it is important to note that the functions hj are nonlinear. So that near the bifurcation, that is for Ae 
small, 

ar, Pr^ 0(\/Ae) but a,; ^ a^o + 0(Ae) and Pi = fiio + 0{Ae) 

In other words it is consistent to consider an approximation only to order V Ae in which case 

Pr = 2pMhAcOs{L0ht + (l)) (39) 

ar — u;Asm{LJht -\- (p) — cos{LUht (j)) (40) 



3 + 2Ae 

" ^/ 8 3« + 47 ^^^^ 



K + 7 2^2k(k + 7)Ae 
4 k(3k + 47) 



(42) 



In figure [T] the actual limit cycles are shown in red and the approximations from center manifold theory in blue 
dashed lines. The limit cycles are shown for various values of the bifurcation parameter Ae and specific values of k 
and 7. 
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QUANTUM NOISE. 



The semiclassical analysis gives us an understanding of the relevant dynamical structures including fixed points 
and limit cycles. However to understand the role of noise in the system it is necessary to proceed to a full quantum 
treatment. While an exact time dependent solution is not available, considerable insight can be gained by considering 
the linearised theory near the fixed points. Such an approach must necessarily fail at the fixed points themselves as 
by definition the linear terms in the dynamical equations vanish at these points, and the eigenvalues of the linearised 
analysis go to zero. Nonetheless such a approach has proved very effective in understanding similar systems in quantum 
optics and we except the same to be true for opto-mechanical systems. 

We follow the treatment of Drummond et al. based on a quantum stochastic differential equations in the 

generalized P-representation. In this approach the quantum operators a, a' , b, are replaced by independent complex 
variables, a,a\ (3, (3'' . Note that the dimension of the phase space has doubled in this representation. This is the 
price one must pay to follow quantum noise for anti-normally ordered moments. The semiclassical equations arise by 
neglecting noise and restricting to the subspace at _> q,* ^ _> ^* _ \Ye will only consider the case of zero temperature 
in this paper, n — 0. 



d_ ( p 

dt \(3^ 

dt [ at 



2ixa 
-2ixa^ 



If we linearise around the semiclassical fixed points (i.e. prior to the Hopf bifurcation e < e/i), the stochastic equations 
of motion are 



d_ 

dt 



5a 
\ 6a^ J 



2 2ixao 2ix(3o 



( 

-2ixal 
2zx/3o 
-2zx/3o* 



\ 



-2ixPo 








which we write as 



dt 



/ 2ixao \ 

-2ixa^ 


\ y 



[Sa] = -A [Sa] + D^/' [qit)] 



2 

1/2 



( 5f3\ 

6(3^ 
Sa 

J \Sa^ J 

/m{t)\ 
vlit) 

\4it)J 



The steady state spectrum for normally ordered moments are defined by 

^'^^^^h l^e— (a.(Oa,(t + r)),_ 



dT 



which is given in terms of the linearised drift (A) and diffusion (D) as 11 1 



S{uj) = ^ {iLvI + A)~^ D {-iljI + A^) \ 



(43) 



(44) 



(45) 



where / is the identity matrix and the superscript T denotes the transpose. 

Direct observations can be made on the output field from the cavity so we look for signatures of the Hopf bifurcation 
in the spectrum of the output field amplitude. In figure [5] we plot the noise spectrum for the optical field amplitude 
for values of mechanical driving approaching the Hopf bifurcation. The spectrum becomes more sharply peaked on a 
frequency uj^ which is equal to the imaginary eigenvalue at the Hopf bifurcation, and given by Eq. (j33[) . After the Hopf 
bifurcation, the real parts of a,f3 are no longer locked to zero but oscillate according to Eq. (36,37). The quantum 
noise will cause phase diffusion around the limit cycle. 

We can estimate this phase diffusion rate using the method of Zhu and Yu[l2] for the response of a van der 
Pol oscillator to white noise in the limit k » 7. Using the linear transformation from a,P to u,v variables to 
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FIG. 2: Linearized spectrum for [S'33(a;)|, for k = 1.0, 7 — 0.1, with three different cases of e = 0.01, e — 0.05 and e = 0.13. 

also transform the noise, we find that both u and v are subject to Gaussian white noise with diffusion constant 
approximately given by s = 1/k close to the Hopf bifurcation. The phase on the limit cycle then undergoes a a 
diffusion process. 



where dW is a Wiener process, and A is the amplitude of the limit cycle given approximately by Eq. (|34p . The phase 
diffusion constant, for k >> 7, is thus 



Like a laser, the phase diffusion rate decreases as the limit cycle becomes bigger. 
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